Evolution to a singular measure and two sums of Lyapunov exponents 
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We consider dissipative dynamical systems represented by a smooth compressible flow in a finite 
domain. The density evolves according to the continuity (Liouville) equation. For a general, non- 
degenerate flow the result of the infinite time evolution of an initially smooth density is a singular 
measure. We give a condition for the non-degeneracy which allows to decide for a given flow whether 
the infinite time limit is singular. The condition uses a Green-Kubo type formula for the space- 
averaged sum of forward and backward-in-time Lyapunov exponents. We discuss how the sums 
determine the fluctuations of the entropy production rate in the SRB state and give examples of 
computation of the sums for certain velocity fields. 
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PACS numbers: 47.10.Fg, 05.45. Df, 47.53.+n 

Dissipative dynamical systems are one of the main 
paradigms of non-equilibrium physics [![. Within the 
paradigm the system is represented by a point in a 
d— dimensional "phase space". The velocity of this point 
is the local value of a prescribed smooth velocity field. 
The latter is compressible so the phase-space volume and 
the Gibbs entropy are not conserved. This setting allows 
to describe effectively non-equilibrium situations with ex- 
change of entropy with the environment. An example of 
such a system, that played an important role in the de- 
velopment of understanding of physics far from equilib- 
rium, is the Lorenz oscillator [2|. This system clarifies 
a rather general feature of the dissipative systems - at 
large times the system's trajectories asymptote a zero- 
volume set in space (attractor). The stationary phase- 
space measure (sometimes given by "the SRB measure" 
[3)@41l) n st is supported on the attractor and is singular: 
it is zero everywhere except on the attractor where it is 
infinite. Here we consider a general dissipative system 
and study the formation of singularities, both infinities 
and zeros, as a result of evolution of an initially smooth 
density. We show that whether infinities or zeros are 
formed is decided by the degeneracy of two time inte- 
grals of a pair-correlation function of velocity divergence. 
The divergence is sometimes referred to as the entropy 
production rate. These integrals correspond to the ex- 
act Green-Kubo type formulas for the sums of forward- 
and backward-in-time Lyapunov exponents respectively, 
where the second formula generalizes the previous result 
in Q . The same quantities are shown to characterize the 
fluctuations of the entropy production rate in the steady 
state. A explicit formula for the sums' difference is pro- 
vided. An essential feature of our treatment is the use of 
Lagrangian trajectories approach that allowed advance- 
ment in a whole range of problems recently, see [7| and 
references therein. The results apply to flows in a finite 
domain that either do not depend on time explicitly or 
depend on time but are stationary with respect to statis- 
tics defined by the spatial averaging (in the latter case 
n st evolves in time and is stationary only statistically). 

Generally a Lyapunov exponent is a non-trivial func- 
tional of the velocity that does not allow for a treatable 
expression in terms of that velocity. As our Green-Kubo 



type formulas show, exceptionally, the sums of Lyapunov 
exponents can be expressed efficiently in terms of the ve- 
locity. The formulas allow a much better control over the 
sum than over any other combination of the Lyapunov 
exponents. Since the sums play a crucial role in determin- 
ing the long-time behavior, then the ability to evaluate 
them is very important. The use of the Green-Kubo type 
formula for the sum of Lyapunov exponents was made in 
the context of studying the behavior of the pollutants 
in the wave turbulence, see d, [Io[ and also [ll|, and in 
analyzing the limit of the weakly compressible dynamics 
Here, as an example, the formulas are applied to a 
random velocity with short correlation time. 

We consider a generic dynamical system governed by a 
deterministic velocity field v(r) defined over some "phase 
space" with a d— dimensional coordinate r. It will be as- 
sumed that the total volume of the available phase space 
is finite and it is preserved by the flow (set equal to unity) . 
The trajectories in the phase space are defined by 



d t q(t,r)=v[q(t,r)}, q(0,r)=r. 



(1) 



It is assumed that v(r) is smooth so that the solution 
to the above system is well-defined and unique. We will 
consider both t > and t < in the above equation. 

We study dissipative systems for which u>(r) = V-v(r) 
is generally non- vanishing. If the trajectories approach a 
zero-volume attractor at large times, then, clearly, under 
the dynamical evolution the volume of finite regions of 
the phase space decays to zero (this does not contradict 
the conservation of the total volume, see the third page). 
In particular, the volume of initially infinitesimal regions 
approaches zero at large times. It is thus instructive to 
consider the evolution of an infinitesimal volume V(t) 
located initially near the point r. This obeys [l3j . 



d\nV 
dt 



u[q(t,r)], In 



V(0) 



Lj[q(t',r)]dt',(2) 



which holds as long as the linear dimensions of the volume 
are much smaller than the scale r\ of smoothness of v. 
The ratio V(t)/V(0) is the jacobian J(t,r) of q(t,r), 



V(0) 



Or 



exp 



u[q(t',r)]dt' 
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It determines the sum of the Lyapunov exponents Ai as 



^Xi(r)= lim - In ^4= lim - 

' V ' t^oo t V(0) t^oo t 



u>[q(t',r)]dt'. 



Thus Ai < would make one expect that initial vol- 
umes decay exponentially to zero at large times. This is 
however beyond this analysis: for any finite initial vol- 
ume, however small, one of the dimensions will eventually 
become larger than r\ and Eq. ([2]) will break down. Here 
we assume positivity of the principal Lyapunov exponent 
Ai (for Ai < the evolution transforms the phase space 
into a collection of point masses) . 

To address the global evolution of the phase space, we 
consider the phase space density n(t,r). Its large-time 
limit at a fixed point r is determined by the limiting value 
of the ratio of the volumes for evolution backward-in-time 
and not forward-in-timc, 



^ XT (r) = lim - In V,. / = - lim 

' » V ' T->ooT V(0) T^oc 



[q(t,r)]dt 
T 



,(3) 



where A~(r) are the Lyapunov exponents of the 
backward-in-time evolution q(—t,r). Density, like any 
transported field Q , is governed by the trajectories prop- 
agation backward (and not forward) in time. To find 
n(t, r) one needs to find the trajectory that arrives at r at 
time t or, equivalently, consider the evolution of the tra- 
jectory backward-in-time starting from that point. The 
Lyapunov exponents for backward-in-time evolution were 
discussed in this context in [8j. We give a self-contained 
presentation avoiding some details not important here. 
The density obeys the continuity (Liouville) equation 



dtn + V • [nv] = 0, 



(4) 



describing mass conservation n[t, q(t, r)]J(t, r) = n(0, r). 
Evolution of density in a moving frame q(t, r) obeys 



lim i ln«^ 

t->oo t n(0, r) 



2 Mr). 



(5) 



To discuss the evolution of density at ia fixed point in 
space we set up the initial condition in the remote past 
at t = —T and study n(r) = n(0, r) at T — > oo, 



n{r) = n[— T, q(—T, r)] exp 



n(r) 



-T 



[q(t',r)]dt' 



lim — In - r 



(6) 
(7) 



Both the logarithmic growth rate of density in the moving 
frame and at a fixed point are determined by the time- 
averages of ui, however the former by the average with 
respect to the evolution forward-in-time and the latter 
by the evolution backward-in-time. These evolutions are 
different and generally ^ A^ ^ ^ A», see [|[ and be- 
low. The difference is made explicit with the help of the 
Green-Kubo type representations for the sums. We use 



[v{r)-V}q(t,r) = v[q(t,r)} 



(8) 



as can be proved by time differentiation [6|. Note that 
if one considers trajectories that issue from the initial 
points r and r+v(r) At then at time t they will be located 
at q(t, r) and q(t+At, r) because the second trajectory is 
the trajectory issuing from r shifted by the time interval 
At. Thus the distance between these trajectories is on 
the one hand dq/dr times the initial separation vector 
v(r)At and on the other hand it is q(t + At, r) — q(t, r) w 
v[q(t,r)]At. Comparing the two expressions one gets 
Eq. ((HJ. Using the latter equation one finds 



uj[q(t,r)]=Vi[q(t,r)] 



dt 



xvj(r) 



duj(x) | 

dxi 



=q{t.r) 



dxi "»=«(*•»•) 
d_ 
drj 



vj(r)- 



d rj 

cj[q(t,r)}. (9) 



We now average over some spatial volume O, 

d f dv f dv 

Jt J^[q{t,r)] = -J-u[rU q{ t,r)] 

+ Ju[q(t,r)]v(r)-dS, (10) 

where the last term is the surface integral over the region 
boundary. Taking to be the whole space the integral 
over the boundary vanishes and one finds 

j t Jdro J [q(t,r)] = -(c(0)u J (t)), (11) 

(w(0)u>{t)) = druj[r}Lj[q{t,r)]. (12) 



The above relation is similar to the relations with 
transient-time correlation functions as described in [l5[ 
and references therein. Integrating Eq. (jll[) we find 



dru}[q(t,r)] 



(w(0)w(t'))df, (13) 



where we used that the LHS vanishes at t = as an 
integral of the derivative. Assuming the integral on the 
RHS converges in the limit \t\ — > oo (remind that t above 
can be both positive and negative) we conclude that, 



//>±oo 
dru[q(t,r)] = - (u(0)u(t)}dt. (14) 
Jo 

Averaging over time, exchanging the order of integrals 
and designating spatial averages by angular brackets, 

J2 x y=~ J o MOM*))*, (is) 

EV>=- f M0)uj(t))dt. (16) 

J — oo 

The above relation depends only on the assumption of 
convergence of the integrals on the RHS. If the conditions 
of the Oseledets theorem [TJ| are satisfied, then X) M r ) 
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and X i (r) are constants (// and /i respectively) for 
almost every r and Eq. (1161) implies 



fi = I (u(0)u(t))dt. 
This is the counterpart of the relation for /i, 



(u(0)ui(t))dt, 



(17) 



(18) 



derived in [6j. The result can be generalized from time- 
independent velocity v(r) to time-dependent velocity 
v(t,r) which is stationary with respect to the statistics 
defined by spatial averaging. One demands equalities like 



J druj(t, r) exp 
J drui(0, r) exp 



- / w[t',q(t'\t,r)]dt' 
o 

o 

u[t',q(t'\r)]dt' 



(19) 



where q(t'\t, r) is the trajectory that passes through r at 
time t. The proof follows the lines of [6| exactly and does 
not need to be repeated here. The formulas (fTT|) - (fT5|) for 
time-independent velocities are a particular case of the 
Kawasaki representation [l5|, [l6[ . 

Based on Eqs. (fT5| - (fT6|) one expects (EAi) and 
(J^A~) are non-positive. However this does not fol- 
low from the non-negativity of the spectrum as the av- 
eraging is not with respect to the invariant measure. 
Also generally (w(0)u;(i)) is not an even function of t 
and ^ A<j ^ \ r • ^ was proved in [f| H, [I?} that 
(E Aj) < 0. Since A~ are the Lyapunov exponents 
of the time-reversed velocity field —v(—t,r), then also 

(E K) ^ °- Note (E K) > would mean E A r (r) > 

for a finite volume, implying by Eq. ([7} infinite mass in 
that region and contradicting mass conservation (thus 
(E A7) — implies ^ A~(r) = almost everywhere and 
the same for E Ai)- 

In the non-degenerate case (Y^, A-j) and (J] A^ ) are neg- 
ative. An important case with degeneracy is where v al- 
lows for a smooth solution to the stationary continuity 
equation V • [n^v] = for the invariant density n st . Here 



u[<l(t,r)] = -wVlnn a t|,(t it .) = -—\nn at [q(t,r)]. 

Thus u>[q(t,r)} is representable as a time-derivative of 
a bounded function. It follows that its time-averages 
EAi(r) and J2^7( r ) vanish (also n st is the same for 
backward- and forward- in-time dynamics [l7|). In the 
generic case of \ t ) < and A~) < Eqs. (fT7|) - (fTH|> 
give complementary descriptions of the formation of sin- 
gularities in the limit of infinite evolution time. The in- 
equality A^ ) < implies the density decays to zero in 
a finite region of space. (Thus roughly (J2 A~) is the rate 
at which trajectories leave the regions outside the attrac- 
tor.) Complementarily, (^2 Aj) < implies the density 
becomes infinite for almost every trajectory starting in a 



finite region of space. (This corresponds to the accumula- 
tion of different trajectories on a zero- volume attractor). 
Nothing seems to forbid (EAi) = at (E A^) < or 
(J2 A~) < at Aj) = (say the rate of density van- 
ishing can be non-exponential). Thus to decide whether 
a given compressible flow will produce singularities out 
of initially smooth density one has to check degeneracy 
of /"jwCOM*))* and / -°>(0Mt)>£fc. 

We analyze the difference between (J2 K) and (J2 A~). 
Using Eq. |3]) we express (^ A^~) via the forward-in-time 
evolution as 

y3 A r) = _ J im / % f dru[q{t,r)} =- lim 

f T dt f \ f T 

' — / dxoj[q(t, x)] exp / oj[q(t ,x)]dt 
o T J J 



(20) 



where we changed variables x = q(—T,r) and shifted 
the time- variable. If uj[q(t, r)] has a finite correlation r c , 
as defined by spatial averaging, then the cumulants of 
oj[q(t', r)]dt' are proportional to T at T > t c and 



7(6) 



lhn (1/T)ln J dr exp 



[q(t',r)]dt' 



is a non-trivial function, as follows from the cumulant 
expansion theorem [18|. From Eq. ([2T)| and the conser- 
vation of the total volume 7(1) = we have 



7 



'(°) = (E A *>< 7(i) = -(E A 



(21) 



Holder inequality implies 7(6) is a convex function which 
together with 7(0) = 7(1) = implies 7'(0) < and 
7'(1) > 0. This is an alternative proof (EA-i) and 
(J2 A~) are non-positive. Further, at t 3> r c the proba- 
bility density P(x, t) of a "sum" x(t) = / * u[g(f, r)]dt' jt 
of a large number ~ t/r c of independent random vari- 
ables has a universal form P{x,t) ~ exp[— tS(x)]. Here 
the "entropy" function S(x) is the Legendre transform of 
7(6). It is positive everywhere except for x = fj,, thus de- 
scribing the exponential decay of the volume fraction of 
regions for which x(t) deviates from its long-time limit /i 
(here constant ^ Aj(r) is assumed). The volume conser- 
vation relation J exp(t[x — S(x)]) ~ 1 is determined by x* 
giving minimum to x — S(x) and obeying X* = S{x*) due 
to 7(1) = 0. While at t — > 00 the volume of the density 
support generally decreases, at any finite t it equals one. 
Evolution of a fraction of space exp[— tS{x*)] expands it 
by the factor exp [fx*] and its volume roughly equals one. 
Often, like for the Galavotti- Cohen relation [l9j . one 

identifies a = — L u>[q(t, x)]dt/T as the entropy produc- 
tion rate in the non-equilibrium steady state described by 
n s t (if it exists). Fluctuations of a in the steady state are 
described by the large deviations function obtainable as 
the Legendre transform of 7(6), 



lib) 



lim — In 

T^oo T 



t (r)drexp 



u[q(t',r)}dt' 



4 



Let us show that the equality between 7(6) and 7(6) 
holds quite generally, cf. @, l20l] . We introduce a time 
i* obeying < f» C T and make a change of variables 
x = q(t*, r) in the definition of 7(6). We find 



7(6)= lim — In / n(t*, x)dx exp 



•T-t, 

/ w[q(t',a;)]di' 



lim —In / n(t*, x)dx exp 6 / cjfqfi', :c)]di' 
T^oo T J Jo 



where n(t, x) obeys Eq. @ with n(0, cc) = 1. Wc shifted 
the time variables in the integral and used i» <C T 1 . Tak- 
ing now the limits t* — >■ 00 and T — > 00 at f» <C T wc 
obtain 7(6) = "/(b). The main assumption is that the 
density n(i*, a;) tends to the stationary measure at large 
t*, cf. [12]. Thus the relations derived previously for 7(6) 
also characterize the fluctuations of the entropy produc- 
tion rate in the steady state. All the above relations can 
be generalized to random, time-dependent velocity 0, Q . 

The difference between (X)Aj) and (J^A~) is due to 
non-even part of (io(0)uj(t)). This vanishes for time- 
reversible statistics. One can derive a series represen- 
tation for the difference in terms of different time corre- 
lations of u) on the trajectory q(t, r). Change of variables 
x = q(—t,r) in the definition ([121 of (cj(0)ui(t)) gives 



<jj(Q)ui(— 1)}= / drbj{r)oj[q(t, r)\ exp 



uj[q(t' ,r)]dt' 



Expanding the exponent into series we find (to = t) 



[E A < - E v] ) = E /°° * / /* B_1 dt n 



drw(r)uj[q(t n , r)]..oj[q(ti,r)]uj[q(t, r)}. 



Eqs. ([L7|) - (fl"8)) allow to compute space-averaged sums of 
the Lyapunov exponents easier than the Lyapunov ex- 
ponents themselves. The latter involve correlation func- 
tions of all orders, not just the second one. For exam- 
ple consider a random velocity field short-correlated in 
time. Then J (uj(0)uj(t))dt is determined by times where 
q(t,r) » r, 



E A ') = (E K ) = "(1/2) / MO, r)a;(t, r))di, 



where we used that the integrand is an even function 
of t. For the so-called Kraichnan model v is a Gaus- 
sian velocity field with zero mean and pair-correlation 
(v a (ti,rx)vp(t2,r 2 )) = [V 5 a p - K a/3 (r 2 - ri)]S(U-ti). 
We find = ^ = (-l/2)\/ a V p K a p(Q), seee. g. @. The 
formulas for other combinations of the Lyapunov expo- 
nents are more cumbersome Q. For isotropic ensemble 
K a/3 (r) = D[(d+1 - 2T)5 af) r i + 2(dT - l)r a rp], where D 
is a constant and < T < 1 is a measure of the flow com- 
pressibility, one obtains /1 = /i~ = —dDT(d — l)(d + 2). 
Another use of Eqs. ((T7| - (|18p is for small-amplitude 
waves such as surface waves in the ocean. Here higher 
powers of v are relatively small. The calculation in [9l-fll| 
showed that up to the fourth order in the wave amplitude 
fi = 0. This is while Ai is of the fourth order in the wave 
amplitude. Thus, though the surface flow is a generic 
flow with compressibility of order one, there is a degen- 
eracy that cancels the long-time effects of compressibility, 

Aj/Ai <C 1. Thus how to explain clustering of pollu- 
tants on oceans' surfaces is unclear. Now, the difference 
between and /i raises the possibility that fi~ does not 
vanish in the fourth order in the wave amplitude. The 
study of this possibility is a subject for future work. 
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